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Abstract 

We reanalyze the recent computation of the amphtude of the Higgs boson decay 
into two photons presented by Gastmans et al. [1,2]. The reasons for which this 
result cannot be the correct one have been discussed in some recent papers. We ad- 
dress here the general issue of the indeterminacy of integrals with four-dimensional 
gauge-breaking regulators and to which extent it might eventually be solved by 
imposing physical constraints. Imposing gauge invariance as the last step upon R^- 
gauge calculations with four-dimensional gauge-breaking regulators, allows indeed 
to recover the well known — )• 77 result. However we show that in the particular 
case of the unitary gauge, the indeterminacy cannot be tackled in this same way. 
The combination of unitary gauge with a cutoff regularization scheme turns out to 
be non-predictive. 
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1 Introduction 



Recently some attention lias been brought back to the VT-loop contribution 
in the calculation of the if — )■ 77 amplitude because of a result presented 
by Gastmans et al. PP] turning out to be at odds with the renowned one 
of Refs. pE]. It goes without saying that, if correct, the result in Refs. [T||2] 
would have had relevant consequences for the ongoing Higgs boson searches 
at the LHC. 

Starting from the observation that the full amplitude if — )■ 77 is free from 
ultraviolet and infrared singularities, Gastmans et al. performed their calcu- 
lation in four dimensions with no regulators and used the unitary gauge to 
consider only the physical degrees of freedom. A gauge invariant amplitude is 
obtained with the 'Dyson subtraction' [5|6] . leading to 
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This amplitude, which happens to vanish in the mw/triH limit (contrary 
to the standard one), would imply a reduction of the decay width r(ii' — )■ 77) 
by about 50% for itlh — 120 GeV, with respect to what found in Refs. [31'^. 

The standard ii 77 amplitude was computed in 't Hooft-Feynman gauge 
with dimensional regularization [3J , with background field methods [1] , and in 
unitary gauge with renormalization group analysis [4J. It reads as 

2 

= (A \! [2 + 3r + 3r(2 - r)/(r)] {k^- k.g^^ - k!^k'()e,{h)e,{h){3) 

Gastmans et al. casted some doubts on the reliability of using dimensional 
regularization and prefer their result adding the motivation that it would 
respect some m\y/mH — )• 'decoupling limit', which has indeed no reason 
to hold true as explained in Refs. |7|8|9|T0] in the framework of equivalence 
theorem [llj . 

The results of Refs. ^Iii2j have been criticized by a number of recent pa- 
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The criticism concerns the absence of regulators, leading to ambiguities in the 
intermediate steps of the calculation, the use of Dyson subtraction and the 
reference to the Appelquist-Carazzone theorem [TB] to justify the decoupling. 

The amplitude if — t- 77 has been calculated in several ways, all confirming 
the old result of Refs. [3|3] , as follows: 

• The authors of Ref. [8] redo the calculation in four dimensions in the unitary 
gauge with a gauge- invariant regularization method (Pauli- Villars like) [IB] . 
They stress the importance of having set an explicit regularization scheme 
to control finite terms which guarantee gauge invariance through all inter- 
mediate steps of the calculation. The authors cross-check the calculation 
with an independent one in dimensional regularization and underscore that 
no renormalization condition should be applied in the presence of only finite 
terms. Thus they conclude that the calculation presented by Gastmans et 
al. must be wrong because it is finite and not gauge-invariant (essentially 
because of the lack of regularization of divergent integrals, finite terms, 
relevant for the gauge invariance of the final result, have been missed). 

• The authors of Ref. perform the calculation in dimensional regulariza- 
tion both in the unitary and in gauges (the same as in Ref. [ID])- It is 
pointed out that without any regulator the coefficient of g^'' , arising upon 
four-dimensional symmetric integration in renormalizable gauges, is an in- 
determinate form of the kind 00 — cxd - also responsible of the breaking 
of gauge invariance - and that in the unitary gauge the same happens to 
the coefficient of /Eg/ci- With reference to Jackiw |T7], the authors stress the 
importance of symmetry requirements (in this case, gauge symmetry) to 
solve the ambiguities in the finite terms occurring in loop calculations. 

• H.S. Shao et al. [12] perform the calculation using a four-dimensional mo- 
mentum cutoff regularization both in 't Hooft-Feynman and unitary gauges. 
Within the latter they obtain the same result of Refs. p|2] starting with a 
particular routing of momenta. They also get the terms to be added after a 
shift in the loop momentum: these contributions sum up to zero when the 
momentum choice of Refs. p!p] is adopted. Performing the calculation in 
the 't Hooft-Feynman gauge, the authors recover the same gauge invariant 
result as the one obtained in dimensional regularization [3] by subtracting 
the countribution of all diagrams (evaluated at ki = k2 = 0) except for the 
(f)~(f)~W~ loop, the one out of the three which is independent of /ci, k2. The 
result is independent of the loop momentum choice because the divergences 
are only logarithmic. 

• F. Bursa et al. [H] perform the calculation of the ii — )■ 77 decay amplitude 
using a (gauge invariant) spacetime lattice regulator and obtain very good 
numerical agreement with the decay amplitude evaluated with dimensional 
regularization. 

Summarizing, in Ref. [S] it is highlighted that the problem of Refs. [T]B] resides 
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in the absence of regulators. In Ref. [12] however it is shown that the use of 
cutoff regularization in unitary gauge leads to confirm the result by Gastmans 
et ai. Thus we might observe that, if there is a problem in the latter calcula- 
tion, it is not in the lack of a regulator but rather in the combination of the 
unitary gauge with the use of non-gauge invariant regulators. As the cutoff 
regularization has been widely used in the literature, we explore further its 
connection with the unitary gauge: in particular we attempt to get a deeper 
understanding of the result presented in Ref. \12\ . 

In Section [2] we focus on the critical integrals at the core of the calculation of 
the — 77 amplitude. We stress, with some elementary examples, that if a 
cutoff regulator is chosen, the definition of the integration boundaries is part 
of the regularization scheme itself. In the paper by Jackiw [ITj it was clearly 
discussed how the critical integrals we have to deal with are indeterminate, as 
long as we do not use a regularization scheme which preserves the full symme- 
try of the theory. Some two-dimensional examples are left in the Appendix |X] 
This would be enough to close the argument here. Yet we show how the inde- 
terminacy in cutoff regularization (Section [3]) can lead to the correct result if 
an appropriate finite renormalization condition is used. 

However in the case of the unitary gauge we will show this is not possible. In 
Section H] we extend the discussion to gauges. Our conclusions are left to 
Section [51 



2 The Critical Integrals in the — t- 77 Amplitude 

At the core of the problem of the i7 — )■ 77 amplitude calculation there is the 
calculation of the integral 



where = ■m^r — XiX2m'jj, and Xi,X2 are Feynman parameters. According 
to Gastmans et ai, performing the integral in four dimensions with symmetric 
boundaries 0, we can substitute l^li, — )• jPg^u-, leading to I^y = 0. Here we can 
appreciate the difference with respect to dimensional regularization (DREG), 



^ Since the integral in Eq. (jlj) does not depend on any external momenta, for 
tensor invariance it must be 7^,^ = Ig^u', by saturating both sides with g^^'^ , we have 
IpLvQ^'^ — ^ 4/, t^g^u — ^ 4/^ and /^/jy — )• then we can solve with respect to /. We 
have the same result if we substitute l^ly — )• ^f'. 
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where it is found 



/£^^^ (-) = f d-l " 3 = -^g,.^- + 0{n-A) (5) 



Gastmans et al. conclude that I^,^^^ (n) must have a discountinuity in 
n = 4, thus mining the foundations of the DREG technique stating that inte- 
grals are not analitic in n dimensions. 

Let us start from the four-dimensional integral in Eq. (j4]). After Wick rota- 
tion and rescaling / — )■ l/M, we get 



V = ^ / dj (6) 



To simplify the discussion, let us focus on the case /x = 1, = 1 



p - ap 

hi = i I dj 7— — ^ =t I dj Fn (/) (7) 



The integrand is not a summable function: it is not positive everywhere in the 
domain of integration, and J d^ \Fu\ = oo, which means that the integral is 
not defined per se - the value depends on how the boundary is chosen to behave 
at infinity. Let us therefore compute the value of the integral over different 
integration domains with different behaviors at infinity: we will observe how 
the integral may assume every finite value, and even diverge. 

As a first example, let us consider a 'spherical cutoff' in the sense described 
below. In polar coordinates, we write 



Iu = i f dl f dn^ (l-A cos^ e) 

Jo (l + Pf J ^ ^ 

= iAn f dl ^ ^ f de sin^ 6 (l - A cos^ e) = 

Jo (1 + Pf Jo ^ ^ 



A is a dimensionsless cutoff, being / a dimensionsless integration variable. 
The angular part vanishes, so there are no problems with the logarithmic 
divergence of the radial part. Actually, every integration domain which has 
the li f-)- ±lj symmetry, leads to an identically vanishing integral. 



^ Which amounts to d^l = id^l, g^y — >• —5fj,u, P — >• P^ and — )■ / 
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h 

Fig. 1. We integrate Fn = — — over an elliptic domain. Because of cylindrical 
symmetry, the graphic is the same independently of Ik = h^h^h- The darker is the 
background, the larger is the Fn value. The boundary is solid when Fu > 0, dashed 
otherwise. Since we broke the spherical symmetry, the ellipse bounds a negative part 
which is larger than the positive one: the integral does not vanish. 

As a second case we choose a non-symmetrical domain of integration. For 
example, let us integrate Fu over the elliptical domain + I2 + I'l + If < A'^ 
(see Fig. ^ 



1 - (4 + 3e) cos^ e 



= " 2^ + «lA^j (^) 

The integral in Eq. ([9]) can assume different finite values as a function of e. 

Choosing asymmetric boundaries, we even loose tensor invariance, obtaining 
a 4 X 4 matrix of unrelated, indeterminate terms. This translates into the fact 
that If^i, is no longer proportional to 6^^, as it should be the case (see Eq (|6])). 
We seek an appropriate choice of the boundaries for all the terms in the I^u 
matrix in such a way to recover a 6f^i, structure. We can therefore compute the 
Ijj,,^ entries by choosing the same asymmetric boundary on all diagonal terms, 
and generic symmetric boundaries for all off-diagonal terms. In this way, all 
diagonal terms will have the same indeterminate value J, whereas off-diagonal 
terms will vanish. We thus obtain I^^^, = 15^1/ , with I being an indeterminate 
(even divergent) constant. 

In the Appendix |A] we give more details on this and we also consider the case 
of Schwinger regularization. We find that, also in this case, the value of the 
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integral changes with the angular dependence of the cutoff function. Schwinger 
regularization turns out to be particularly instructive as the indeterminacy of 
the integrals is not related to the shape of the integration domain but to the 
behavior at infinity of the integrand function. This is a signal of the fact that, 
in order to tackle the indeterminacy of the critical integrals in this calculations, 
one exclusively needs to set physical constraints as there is no mathematical 
prescription which can univocally determine them. 

Gastmans et al. rely on their finite (equal to zero) result for the integral of 
the type of Eq. (jl]), which follows from a particular choice of the integration 
domain. This also explains why the calculation of Ref . [12] , in which a spherical 
cutoff is explicitely used, leads to the same result found by Gastmans et al. in 
unitary gauge: the choice of the integration domain in Ref. [12] is the same as 
the one implicitely taken in Refs. [TP] . 

The authors of Ref. [9] , on the other hand, underscore the fact that the inte- 
gral (jl]) is an indeterminate oo — oo form: it must be treated with some reg- 
ularization scheme. Gauge-invariance in the final result can be implemented 
either a priori by choosing gauge-invariant regulators (like Pauli-Villars or 
DREG), or a posteriori by applying an appropriate subtraction. We remark 
that in the latter case integrals are not well defined, and their values must be 
considered indeterminate. 

If we choose a sharp spherical cutoff we still get the Gastmans et al. result 
in unitary gauge as was first shown in Ref. [12]. On the other hand, if we 
use renormalizable gauges in a cutoff scheme, we recover the standard result 
('t Hooft-Feynman gauge is used in Ref. [12]). 

Does the fact that two different results are obtained using two different gauge 
choices mean that a cutoff scheme is not to be pursued at all? In the follow- 
ing we mean to show that the problem does not reside in the use of a cutoff 
scheme by itself but rather in figuring out that different ways (spherical, ellip- 
tical etc.) of implementing a cutoff scheme amount to different values of the 
integrals, i.e., to indeterminate coefficients. We can actually use some cutoff 
scheme provided that there is a clear recipe on how to absorb the indetermi- 
nate coefficients arising in the calculation, also restoring gauge invariance. We 
will show that such a recipe cannot he found in unitary gauge. 

After this discussion one might therefore ask why cutoff renormalization works 
in several cases independently of the issue of the indeterminate constants 
discussed above. 

For example, let us discuss the cutoff renormalization of vacuum polarization 
in QED. Since we are not protected by Ward identity, the integral presents a 
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quadratic divergence 



/ 



nf (k) = -Uel 



(10) 



The integral is divergent and sign-undefined; it is ill-defined as the one in 
Eq. (jl]). We could repeat the above considerations to show that we can lower 
the degree of divergence with an appropriate choice of boundary. Anyway we 
simply add a photon mass counterterm and impose that the quadratic term 
vanishes (i.e. that the photon is massless), in order to recover a good gauge in- 
variant theory. The presence of the counterterm assures that the quadratic di- 
vergence disappears whatever value the integral has. In this case, we can forget 
about the indeterminacy itself, and calculate the integral with the boundary 
we prefer. On the other hand, we cannot choose a boundary and claim that 
an integral like the one in Eq. has a finite value and needs no counterterm 
(see Gastmans et al); we must consider expressly the indeterminacy and add 
a counterterm to absorb it. 



3 Considerations on power counting and gauge invariance 

Jackiw [iZl showed that indeterminacy can arise if we use regulators which 
have less symmetry than the theory, and it is not necessarily resolved when 
we restore the symmetry at the end of calculation. On the other hand, if the 
regulator maintains the symmetry, the result will be univocal. We can simply 
understand the case of gauge symmetry: gauge invariant regulators decrease 
the degree of divergence of the integral, making it finite and regulator indepen- 
dent. Indeed, by naive power counting, we know that the amplitude if — )■ 77 
is logarithmically divergent in renormalizable gauges. In gauge invariant regu- 
larizations, we can group two momentum powers in the numerator to extract 
the gauge invariant factor ki - k2g^'^ — kl^ki, so that the amplitude becomes 
finite. In unitary gauge too, we expect the same finite amplitude after a not- 
straightforward cancellation of higher divergent terms. However, this cannot 
be done in cutoff regularization where the Ward identity and gauge invariance 
are spoiled by the breaking of shift invariance. In the latter case the integral re- 
mains indeterminate or at worse logarithmically divergent. The expectations 
by Gastmans et al. to get a finite amplitude which needs no regulator are 
disappointed by the choice of four-dimensional symmetric integration, which 
implicitly uses a spherical cutoff scheme, leading to the breaking of gauge 
invariance and to a divergent amplitude. 

To subtract the divergence in the cutoff regularization scheme and, in gen- 
eral, all cutoff-dependent terms, we need some counterterms. Breaking gauge 
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symmetry, we have the most general lagrangian with all possible combina- 
tions of bare fields and bare couplings, even an ad hoc counterterm of the 
form diriQA h^A^. This is what Dyson subtraction means: hide all divergent, 
cutoff-dependent, non-gauge-invariant terms into a counterterm, which would 
vanish in a gauge invariant regularization scheme. 

We have computed if — >■ 77 amplitude in 't Hooft-Feynman gauge without 
calculating divergent integrals; we find: 



-ki^kl (2 + 3r + 3r (2 - r) / (r)) 



2m]j 



dxi 



1 — Zl 



dxo 



dH 
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^ (/2 - 1 + AX1X2T + ie) 



+ >k^Ml + ^^ + 3r(2-r)/(r) 



(11) 



We remark that the first term (proportional to kl^ki) contains only well-defined 
finite integrals; the second term is indeterminate (vanishing according to sym- 
metric integration as in Refs. [TjE]). With the use of DREG, the second term 
would give \m'\jg^'' {l + , leading to the well known gauge-invariant ex- 
pression. However, let us stay in the framework of gauge-breaking regulariza- 
tions. 

In Ref. [12] a modified version of Dyson subtraction is performed to recover 
gauge invariance. One might even wonder whether Dyson subtraction is al- 
lowed without divergent terms p]. As we have just shown, the integral in the 
second term is probably divergent and in any case cutoff-dependent, so we are 
allowed to add a counterterm and impose gauge invariance as a renormaliza- 
tion condition. In so doing we get the correct amplitude in Eq. (j3]). We would 
have the same expression by using symmetric integration: every value of the 
integral disappears into the counterterm. We are therefore lead to observe that 
the arbitrariness related to the choice of the boundary (or, in general, of the 
regulator) is solved by imposing gauge invariance. 

Why Gastmans et al.^s amplitude is different from the standard one? As shown 
in Ref. [2], in unitary gauge we have another divergent integral 



rl-xi 




L '''' J 


[ dH 







{k^kl - ki ■ kiQ^") f - 2k\ {k2 ■ I) l^" 

~2kl^ {ki -1)1" + 2 {ki ■ k2) + 2c/^'^ {ki ■ I) {k2 ■ I) 
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1 



{P - myy + XiX2mjj + ie) 



(12) 



By symmetric integration the integral vanishes, whereas dimensional regular- 
ization leads to ^']3j^jt;q ('^) = ^^r^ {k2k^ — ki ■ k2g^") + {n — 4). The integral 
has the same indeterminate behavior of the former one: the value depends on 
the choice of the boundary. We can say that A = Jkl^k'^ + J'g^"^, with J 
and J' indeterminate constants. While in the Eq. ffTTl) the tensor kl^k'^ has a 
well-defined finite coefficient, and we can tune the rest of the amplitude on it, 
in unitary gauge this coefficient is indeterminate, possibly divergent: we must 
then add another counterterm SgoAho (S'^Aq)^ to absorb the divergence. 

We have now two counterterms and we need two renormalization conditions 
to fix the arbitrariness. The only Dyson subtraction (which means imposing 
gauge invariance) is not enough anymore. The result in Eq. ([1]) is still arbitrary, 
and allows the addition of whatever gauge invariant — ki- k2 g^^ term. 
The other condition could be, for example, the requirement of the validity of 
the equivalence theorem [11] in the limit my/ — )■ 0, or the invariance of the 
amplitude in both 't Hooft-Feynman and unitary gauge: both conditions fix 
the value of the amplitude in Eq. ([1]) to the standard result in Eq. ([3]). The 
indeterminacy is resolved; we recover also the independence of the amplitude 
on gauge choice and regulator choice. 



4 Renormalizable gauges 



The main difference between renormalizable gauges and the unitary gauge 
is in the ultraviolet behaviour of the propagator 



lA^^ = , ^ , (g'-'' - (1 - 2 2 ] 
-m^ \ q^ - ^mfy J 

- 'if q'^q^ 1 

= - T + ^ 2 , 2 (13) 

For finite ^, the propagator is 0(g~^), and leads to a logarithmic divergent 
integral, the same as in 't Hooft-Feynman gauge; on the other hand, for ^ = oo 
(unitary gauge) the propagator is O (g°), leading to a highly- divergent ampli- 
tude. If we use a gauge-invariant regulator, we could simply handle all high 
divergences and show that all ^-dependent terms in the amplitude vanish, 
leading to the same result as in the unitary gauge [9]. In other words, since 
gauge-invariant regularizations give meaning to the integrals, we can carry the 
limit = oo under the integral. 
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On the contrary, in gauge-breaking regularization, divergent integrals have 
ambiguities, and we have to handle limits with care. Let us study the divergent 
behaviour in gauges. The integrals remain logarithmic divergent, so we can 
perform usual Feynman parametrization and shift the integration variables. 
We obtain the expression 



I (27r) mw ^ ^ Jo i=i \ i=i J 

dH ( 4Fr - fg^"") + finite integrals (14) 

{12 _ M2 + lef ^ ^ ^ ^ ^ ' 

with = (1 + — 1) (x4 + X5 + xq)) — m^j (xs + xq) (xi + X4). Despite 
of the complication of Feynman parameters, we have the same form as in 
Eq. ( ITTl) : the integral can be treated as an indeterminate tensor Ig^'^, while 
the coefficient of ^2^1 is finite and uniquely determined. As in the previous 
section, by Dyson subtraction we can tune the indeterminate constant and 
recover the correct gauge invariant amplitude in Eq. ([3]). 

In the unitary gauge, the behaviour is completely different: the divergent inte- 
gral in Eq. (fH|) should vanish in the ^ — )■ 00 limit, so the divergences in unitary 
gauge must arise from a combination of highly- divergent terms. Indeed, let us 
consider again the calculation of Ref . [2] . Putting external momenta on-shell, 
Gastmans et al. show that high divergences reduce to O (fc^) terms, and then 
get simplified to the logarithmic divergent term in Eq. (1121) and in a quadrat- 
ically divergent term which, with an appropriate choice of loop momentum, 
vanishes for tensor invariance [ll- 

To summarize, the indeterminate behaviour of k2k\ in the unitary gauge arises 
from highly divergent terms which do not appear in gauges. The limit 
— 7- 00 cannot be taken under the integral because integrals are ill-defined, 
and integrand functions are not measurable. The amplitude is arbitrary only in 
the unitary gauge, and divergences apparently disappear only with a particular 
choice of loop momentum. 



5 Conclusions 



We have analyzed the computation of the amplitude if —t- 77 by Gastmans 
et al. [1|2], to understand why it turns out to be different from the standard 

^ An integral of the form f dH ^ — . , . , must vanish because we have no 1- 

° J even polinomial m t 

index constant tensor it can be proportional to. Obviously, it would vanish also in 
symmetric integration. 
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result in Eq. ([3]). Integrals of the form 



/ 




(15) 



are not well defined. We have provided some explicit examples within cutoff 
regularization, obtaining different values by varying integration boundaries. 

In 't Hooft-Feynman gauge and in a cutoff regularization scheme, see Eq. (|TT1) . 
we obtain 



MIL, = '/^ [-k^K (2 + 3r + 3r (2 - r) / (r)) + Ig^ (16) 



where / is a constant which depends on the boundary shape. This makes 
it indeterminate as there is no physical prescription on the choice of the in- 
tegration boundary shape. On the other hand, the use of a gauge invariant 
regularization scheme provides automatically the recipe on how to evaluate 
the integrals. 

Since the term kf^k'^, in Eq. f|T5|) . has only one finite unambiguos coefficient, 
we are able to solve the indeterminacy by imposing gauge invariance at the 
end of the calculation. 

However we have shown that, in the unitary gauge, both the coefficients kl^k'^ 
and g'^'^ are indeterminate in the sense of integration boundary shape de- 
pendency. Imposing only one renormalization condition (like imposing gauge 
invariance by Dyson subtraction) is not enough anymore. Given the equiva- 
lence of gauges with unitary gauge as ^ ^ oo, the problem we discuss is 
likely related to the exchange of this limit with integral sign for not Riemann- 
summable functions: the coefficient of kl^k, arises from highly divergent terms 
which do not appear at finite values of ^. 

Gastmans et a/.'s expression in Eq. ([T]) is still ambiguous upon Dyson subtrac- 
tion, and allows the addition of whatever term of the form kl^k'^ — k, - k2g^^. 
This arbitrariness can be fixed by requiring the validity of the equivalence 
theorem, or by imposing the equality of amplitudes in unitary and 't Hooft- 
Feynman gauges. In other words, we are able to add terms to Eq. ([1]) in order 
to match the standard result Eq. ([3]). 

The combination of unitary gauge with a cutoff regularization scheme simply 
turns out to be non-predictive. 



,2 
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A Appendix 



For the sake of simplicity in the following we consider a two-dimensional ver- 
sion of If^i, in Eq. 

V = ^ / d2l ^^"^^ ~ ^V'^ (A.l) 



This version of I^u has the same properties of its four-dimensional counter- 
part, namely: i) the integral is superficially divergent as a logarithm, ii) it is 
identically zero for symmetric integration domains. Hi) the integrand function 
has no definite sign. The conclusions we will draw from the following calcu- 
lations in two dimensions remain unaltered in four dimensions: here we just 
avoided superfluous technical complications. 

As we did before, let us start by computing the In term. We realize that In 
can be mapped into an entry of the /12 kind upon a rotation by 45" of I1I2 
axes. Ill = /12 only if the integration domain is rotated accordingly. Since the 
calculations turn out to be simpler using the {12} entry, we will make our 
observations on this case only 

/12 = ^ y ^2^77^^ =1 j d^i (A.2) 



At any rate we remark that the domains of integrations will be chosen in 
such a way that eventually all off-diagonal I^y entries will vanish as to recover 
eventually the 5^i, tensor structure. 



The integrand in Eq. flA.2|) is negative when I1I2 > (I and III quadrant), and 



positive otherwise. In the former case, we bound a domain with two quarters 
of a circumference of radius A; in the latter case we use a square with edge A 
(Fig.lAll). We have 



Ii2 = —Ai I dl TT / c?^ sin ^ cos ^ + 4i / dUdU n 

io (1 + Jo Ao,A]x[0 A] (1 + /f + llf 

Again we get a finite non-zero value. The leading divergences are the same 
in each quadrant, whereas the finite part is boundary- dependent, so that the 
sum does not vanish. 
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Fig. A.l. We integrate Fi2 = ~ ,oh over a mixed boundary. We bound the domain 

(i+t ) 

with a circle when F12 < 0, and with a square otherwise. Since we broke spherical 
symmetry, the integral does not vanish. 



More in general, we can slice in a countable set of bounded regions, in order 
to reduce the integral over the whole to a countable sum of finite integrals, 
i.e., to a series. We can thus use the Riemann rearrangement theorem [18] to 
obtain whatever finite value or logarithmic divergence. 

For example, let us consider all the concentric circumferences with integer 
radius thus slicing into annuli: the integral of F12 over each annulus vanishes 
by circular symmetry. Therefore we slice each annulus into a positive region 
Pk where F12 > 0, and a negative region A^^^ where F12 < (Fig. IA.2I) . We 
therefore have 



r /"""/S fk+l 

Pk= d2l F12 = 4: / sin 6' COS / dl 

J Pk Jo Jk 

1 1 fc2 + 2fc + 2 

~ P + 2A; + 2 ~ P + 1 ^ °^ + 1 



rik = / d2l F12 = -Pk (A.4) 

■JNk 

The Pk form a bounded sequence of positive terms converging to 0. We can 
find that the greatest term of the sequence is pi = M ^ 0.62. Specularly, 
the nfc form a sequence of negative terms converging to zero, bounded by 
rii = —M ^ —0.62. If we the the union of all Pk and Nk, we recover the whole 
M^, therefore if we sum all pk and we recover the whole integral. Since 
J2kPk and Y.k'^k both diverge separately, we must specify the correct ordering 
of terms. We start by adding the first positive terms Pk until we exceed 1 + M, 
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(a) (b) (c) 

Fig. A. 2. Riemann rearrangement. Light gray regions have positive integral pk, dark 
gray regions have negative integrals n^. The absolute value of each region is bounded 
by M ~ 0.62. We can obtain a divergent sum following this algorithm: (a) we sum 
first positive terms pq + . . . + pN^, until we exceed 1 + M, then we can subtract 
no still exceeding 1; (b) we continue adding positive terms until we exceed 2 + M, 
then we can subtract ni still exceeding 2; (c) and so on. We see that the negative 
region becomes smaller and smaller than the positive region, so that cannot cancel 
the logarithmic divergence. 

and then add the first negative term no- Since all > > — M, we have still 
Po +Pi + • • • +PAri - l^ol > 1 (A. 5) 



We can continue adding positive terms until we exceed 2 + M, and then add 
ni, and so on. The resulting sum covers all Pk and regions. The series 
diverges, and so does the integral. 

One might wonder whether we have the same behavior with a smooth cutoff. 
We calculate Eq. (]A.2I) with a Schwinger regulator [19] 



7^ 

2tt 



= -iT(o,^^j J rf^sin^cos0 = O (A.6) 

where T{a, h) is the incomplete Gamma function [20]. Again, the angular part 
of the integral vanishes, and so we do not care about the logarithmic diver- 
gence in the radial part. However if we deform the cutoff giving an angular 
dependency to it, e.g. A — )• A (0) = Aexp (e sin 26*) (Fig. IA.3|) . we obtain 



Ii2 = i I (—2/1/2) I dsse 
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Fig. A. 3. We integrate Fyi = ^'i^A regulating the function with a smooth cutoff 

(i+t ) 

A. In the graphic we show a deformed cutoff A (^) = Ae*^*^™^^. The curve is fihed 
in the region where Fyi > 0, dashed otherwise. Since we broke spherical symmetry, 
the integral does not vanish. 



dB sin^cos^r 



' A2(^) 



-me + O 



1 
A2 



(A.7) 



Also in this case, the value of the integral depends on the shape of the cutoff 
function, no matter if smooth or sharp. 
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